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NACA ACR No. L5B01 . 

NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

ADVANCE CONFIDENTIAL REPORT 

DETERMINATION OF THE EFFECT OF HORIZONTAL-TAIL 
FLEXIBILITY ON LONGITUDINAL 
CONTROL CHARACTERISTICS 
By S. M. Harmon 

SUMMARY 


An Iteration method is given for determining the 
longitudinal control characteristics of a flexible 
horizontal tail. The method permits factors such as the 
actual spanwise variation of elasticity and the aero- 
dynamic Induction effects due to three-dimensional flow 
to be accounted for to any degree of accuracy appropriate 
to a particular case. 

An analysis Is Included of the effects of horizontal- 
tail flexibility on the tail effectiveness, the hinge- 
moment characteristics, and the control-force gradients 
in a dive »»covery for two modern fighter airplanes. The 
effects of variations In speed, altitude, elevator 
stiffness, and center-of-gravi ty movements are considered. 
The results of these calculations for speeds below that 
at which critical compressibility effects occur Indicate 
for the two airplanes significant effects due to the tail 
flexibility. It appears that the location of the flexural 
axis of the stabilizer too far behind the aerodynamic 
center of the tail may cause excessive control forces in 
a dive recovery at high speeds. 


INTRODUCTION 


The design of tail structures for high-speed flight 
requires special consideration of the factors that 
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provide sufficient rigidity In torsion In order to 
ensure satisfactory control and maneuverability for the 
complete speed range. Reference 1 presents an analytical 
treatment of the effect of horizontal- tail flexibility 
on longitudinal control characteristics. The analysis 
of reference 1 is based essentially on the assumption of 
a semirigid tail structure having a linear spanwise 
twist distribution and on two-dimensional section force 
theory. The assumption of a semirigid tail, however, 
does not provide for the establishment of the required 
equilibrium between the aerodynamic and elastic forces 
at every section; consequently there is, in general, no 
assurance of the extent to which the arbitrarily chosen 
twist distribution represents the distortion of the 
actual flexible tail. In addition, the low aspect 
ratios commonly employed on tails produce sigrfificant 
induced effocts on the aerodynamic forces. It appears, 
therefore, that more reliable predictions of the control 
characteristics of a flexible tail could be obtained by 
taking account of the actual spanwise variation of 
elasticity and of the aerodynamic induction effects. 

The present paper presents a method for determining 
the control characteristics of a flexible tail that 
takes account of factors, such as the actual spanwise 
variation of elasticity and the aerodynamic Induced 
effects, to a degree of accuracy appropriate to any 
particular case. The method is based on an iteration 
procedure in which the effect of the tail flexibility Is 
obtained by means of a series formed by the addition of 
the Incremental effects resulting from each Iteration. 

The rapidity of the convergence of this series depends on 
the degree of rigidity of the tall, and the Increments 
for the higher-order Iterations can usually be estimated 
from a knowledge of the values obtained from the preceding 
Iterations . 

In order to Illustrate the iteration procedure and 
to indicate the magnitude of the effects of tail 
flexibility in some tyolcal cases, the present investi- 
gation Includes an analysis for two modern fighter air- 
planes of the effect of horizontal-tail flexibility on 
the tail effectiveness, on the hinge-moment charac- 
teristics, and on the control-force gradients required 
In recovery from a dive. The results of these compu- 
tations are given for sea level and for an altitude 
of 30,000 feet for a speed range corresponding to 
Mach numbers ranging from 0 to 0.72. 
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SYMBOLS 


M pitching-moment contribution of tail about 

center of gravity of airplane, positive when 
airplane noses up, foot-pounds; Mach number 
when used to account for compressibility 
effects 

Ifj tail length, measured from center of gravity of 

airplane to elastic axis of tail, feet 

e distance from aerodynamic center to flexural 

center at a section for the tail, positive 
when aerodynamic center Is ahead of flexural 
center, feet 

p air density, slugs per cubic foot 

V true airspeed, miles per hour 

q dynamic pressure, pounds per square foot 

jl 1.1-67 ) 2 fv 2 ] 

T total torque of tail, positive when stabilizer 

leading edge tends to nose upward (the 
dT 

derivative — represents the torque per 
dii 

unit span at a tail section), foot-pounds 
c wing chord, feet 

b span (of wing, unless otherwise indicated), feet 

S area (of wing unless otherwise indicated), square 

feet 

c e root-mean- square elevator chord, measured behind 

hinge line, feet 

c w mean aerodynamic chord of wing, feet 

y coordinate indicating fixed position along span 

from center line 

r] coordinate indicating variable position along 

span from center line 
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aspect ratio 

fictitious aspect ratio employed in corrections 
for compressibility effects (a V l-M^ ) 

ratio of semiperimeter of ellipse to span of 

airfoil surface, primed to indicate fictitious 
plan form employed in corrections for 
compressibility effects 

two-dimensional lift-curve slope for tail 


section lift coefficient for tail; primed to 
refer to section lift coefficient of a 
fictitious plan form employed in corrections 
for compressibility effects 

geometric angle of attack of the tail, measured 
from zero-lift line at section for assumed 
rigid tail 

angular deflection of stabilizer due to tail 

flexibility, positive when leading edge moves 
upward, degrees 

geometric angle of attack of tail, measured from 
zero-lift line at section for flexible tail, 
degrees (at R + 0^ 

elevator deflection at section for assumed rigid 
tail, positive when trailing edge moves 
downward, degrees 

angular deflection of elevator section due to 
elevator flexibility, positive when trailing 
edge moves downward, degrees 

elevator deflection at section in flexible tail, 
degrees (6 r +0-0) 

change in 6p per unit change in normal 

acceleration in recovery from dive, degrees 
per g 


change In per unit change In normal 

R 

acceleration In recovery from dive, 
degrees per g 

change in elevator oontrol force per unit 
change In normal acceleration, pounds 
per g 

acceleration of gravity, 32.2 feet per second^ 


rate of change of section angle of attack with 
elevator deflection for constant lift at 
section for assumed rigid tail 


rate of change of section hinge-moment coef- 
ficient with section lift coefficient for 
constant elevator deflection for assumed 
rigid tail 


rate of change of section hinge-moment coef 
ficient with section elevator deflection 
of assumed rigid tail in degrees for 
constant section lift 


rate of change of section pitching-moment 

coefficient with section elevator deflection 
of assumed rigid tail in degrees for constant 
section lift 


airplane lift coefficient 

three-dimensional slope of lift curve for 
airplane 


rate of change of elevator deflection with 
airplane lift coefficient for trim for 
assumed rigid tail, degrees 
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rate of change of downwash angle at tall with 
wing angle of attack, degrees per degree 


r— ■- comDressibllity correction factor, where 
VL - M 2 M = Mach number 


compressibility correction factor 


/ EA + 2 \ 

V^E'A* + Zj 


elevator gearing ratio, as obtained with no 
load on tail, radians per foot 

airplane gross weight, pounds 

total hinge moment on elevator, positive when 
leading edge tends to move upward, foot- 
pounds 

... . / H \ 


elevator hinge-mom9nt coefficient 


?°e 2b e 


rate of change of with 6 R as obtained 

for given movement of elevator control 
sti ck 

rate of change of Ch with over tail 

pitching-moment coefficient due to tail about 
center of gravity of airplane (M/qSc w ^ 

rate of change of Cm with 6 R , as obtained 

for given movement of elevator control 
stick 


rate of change of with a t over tail 


rate of change of 6 R for given movement .of 
elevator control stick with a+- u over tail. 


with Cm constant 


f-6C m /da t 
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t(y) 


h(y) 


CTRg^) 


^TRe^) 


dg 

dH 


total' torque transmitted by. the stabilizer 
section at station y, foot-pounds 

/ f*> t/2 dt 


cr »*") 


total hinge moment transmitted by elevator 
section at station y, foot-pounds 
' K t /2 

dh , 

— dt! 
d-n 

'y 


coefficient of torsional rigidity for 
stabilizer at station -q, equal to 


t (ti ) 


de/d-q* 

where d0/dn Is the slope of the 0 curve 
at station -q, pound-feet^ per degree 

coefficient of torsional rigidity for elevator 
at station q, equal to where dg/c 

1 s slope of g curve at station q, 
pound-feet2 per degree 


rate of change of section elevator twist with 
total hinge moment on a loaded half of 
elevator surface a3 measured in static 
tests, degrees per foot-pound 


Subscripts : 
t tall 

w wing 

e elevator 

s stabilizer 

R refers to assumed rigid tail 

0,1, 2, etc. numerical subscripts used to Indicate the 

order of twist iteration 
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PRESENTATION OP METHOD 
Development of Formulas 


The pitching moment due to the tail about the 
center of gravity of an airplane, considered positive in 
the nose-up condition, is given by the equation 


M 




(i) 


where T is the tail pitching moment about the flexural 
axis of the tail, assumed to be positive when the 
leading edge tends to move upward. In conventional 
airplanes, the value of T usually increases the 
elevator effectiveness numerically by about 5 percent. 

If the lifting-line theory of reference 2 is followed, 
the lift coefficient c in equation (1) is givsn as 

a function of the 3panwise coordinate y in the form 


oi t (y) 


= a r 





dc 


l t c t 


djL 


dp 


y - h 


U) 


in which, for the flexible tall, 

a t = a t R + 9 
5 = + $ - 9 

The integral expression in equation (2) represents the 
induced uownwash angle. The determination of the lift 
distribution by means of the lifting-line theory for an 
arbitrary angle-of -attack and chord distribution ha3 
received much attention, and numerous methods 
(references 2, *, and 4) are available for obtaining the 
solution of equation (2) when the functions and 6 

are given. 

The basic consideration in the determination of the 
spanwise twist distributions for the stabilizer and 
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elevator, 0(y) and #(y), for use in equation (2) is 
the establishment, at every section, of equilibrium 
between the aerodynamic and the=elastlc forces acting, on 
the tail structure. In the present analysis, 0 and 
are determined on the basis of the theory of pure torsion 
of tubes (references 5 and £)• Other considerations 
relating to the torsional effects of the axial stresses 
induced by the restraints to the free warping of the 
sections of the stabilizer and elevator and to the effects 
of the bending of the ribs can be accounted for by 
employing the proper parameters and following the 
procedure of successive approximations or iterations 
described herein under ‘'Iteration Method. “ 


The aerodynamic twisting moment for a symmetrical 
airfoil section results from the lift distribution 
contributed by the angle of attack, which acts at the 
aerodynamic center of the section, and the lift 
distribution contributed by the elevator deflection, 
which acts at Its center of pressure. If the section Is 
unsymmetrical, the lift distribution due to camber 
contributes a further Increment to the twisting moment. 
On the basis of the foregoing assumptions, If a 
symmetrical section Is used, the applied twisting moment 
across a section dp of the tail Is 


dT(ri) 



+ 



dp 


( 3 ) 


In order to obtain the torsional moment on the 
stabilizer, the moment acting about the elevator hinges 
should be deducted from the total twisting moment on the 
tall, because the elevator hinge moment is normally 
transmitted to the fuselage through the torque tube. 

The applied twisting moment across a section dp of the 
stabilizer, therefore. Is given by 

dt (p) = dT(p) - dh(p ) (If.) 

where dh(p) Is the elevator hinge moment at the 
section p of width dp and 



dh(p) = 


( 5 ) 
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The expression for oi t in equations (3) and (5) is 
given in equation ( 2 ). 

The total torque transmitted by a stabilizer 
section y is 


t(y) 



( 6 ) 


Division of equation (lj.) by dr) and substitution of the 
value obtained for dt/drj in equation ( 6 ) gives 



(7) 


dT dh 

where — (r|) and — (q) are given in equations i' 3 ) 
dn d-q 

and ( 5 ), respectively. Similarly, the total elevator 
hinge moment transmitted by a section y is 


h(y) 



( 8 ) 


If the boundary condition that the twist is zero at 
the root is assumed, the angles of twist for the 
stabilizer and elevator can be expressed in the form 


e (y) = 


(9) 

#(y) = 


(10) 


The torsional-rigidity coefficients for the 
stabilizer and elevator, respectively, are defined as 
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c T R a (-n) 


t(n) 

dO/dp 


C TR 0 ^^ 


h(-n) 

d#/dp 


where t(ri) and h(p) refer to the total torque 
transmitted by a section p. Substitution for dG/dp 
and d#/d p in equations (9) and (10) results in 


e(y) 



( 11 ) 


#(y) 



dp 


( 12 ) 


where t(p) and h(p), corresponding to t(y) and h(y), 
are given by equations (7) and (8), respectively. 


Equations (2), (7), (8), (11), and (12) express, 
within the limitations of the theory involved, the 
equilibrium conditions at each section between the aero- 
dynamic and elastic torques. iiThen the aerodynamic, 
geometric, and structural parameters expressed in these 
equations have been determined, the three unknown 
variables 0, and c i remain. The simultaneous 

integral equations resulting from the required equilibrium 
condition generally involve complicated functions for the 
three unknown variables. In practical cases, however, it 
has been found convenient to determine the characteristics 
of the flexible tail by working with the integral 
equations through a procedure of successive approximations 
based on an Iteration procedure. 


Iteration Method 

The first approximation to the tail configuration 
is taken as the one corresponding to an assumed rigid 
tall; that is, 0 and $ are both zero, and the elevator 
deflection 6 and geometric angle of attack of the tail o^. 
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at each section are equal, respectively, to 6 r and a tj^* 

Corresponding values of c ltQ (y) are "then determined 

from equation (2). Substitution of these values of cj tQ (y) 

in equations (3) and (5) gives ^(p) and ~(p)» and 

dr) dr) 

the functions tg(y), ho(y), 0i(yK and ^(y) can be 

determined, respectively, from equations (7), (8), (11), 
and (12). The values thus determined for 9^ and fa 

can then be employed in turn to determine, successively, 
new increments in the c j^., 0, and 0 distributions. 

The series resulting from the addition of the successive 
increments permit the determination of the control charac- 
teristics for the flexible tail. 


The general procedure for determining th9 elevator 
effectiveness of the flexible tail is as follows: assume, 

as a first approximation, that the geometric angle of 
attack a-fc(y) is equal to zero and the elevator 

deflection 5(y) is equal to 6r. Compute ci to (y) 

from equation (2). Obtain values for dT/dt) and dh/dr) 
at several spanwise stations from equations (3) and (5), 
respectively, with 6(y) = 6r and cj t (y) = cj^fy). 

Integrate equations (7) and (8) to obtain, respectively, 
t 0 (y) and ho(y). Substitute the values of tQ(y) and 
hg(y) corresponding to tg (p) and hg(p) into 
equations (11) and (12), respectively, and obtain ©^(y) 


and fa(y) as a first approximation to the twist 
distributions. For the second iteration (first twist 
iteration), assume that 6(y) = 0±(y) - 6]_(y) and 


that at(y) = 0i(y) 


and compute the corresponding 


c 


U 


1 


distribution from equation (2). The substitutions and 
integrations in equations (3), (5), (7), (8), (11), 
and (12) with 6 = fa - 0 i and c i t ± cj t in the 


manner described for the previous iteration then provide 
the second twist increments 02(y) and 02(7)* For the 
next iteration, assume 6(y) = fai'j) - 02(y) and 


a-fc(y) = 02(7) and obtain, as described previously, the 
third twist increments, 0^(y) and fa(y) . This 


iteration procedure is continued until the increments 
for cj t , 0, and ft become negligible. 
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The foregoing procedure for obtaining the distri- 
butions .. c.j.t., 0*. and is summarized in the following 

table, which gives the variables employed in each " " 
iterations 


’'■“'■'■'■^^Order of twist 
— -—^iteration 
Variable ^ — — 

0 

1 

2 

5 

0 

0 

01 

e 2 

e 3 

0 

0 

01 

02 

03 

“t 

0 

01 

02 

0 3 

S 

6R 

(01 - 0i) 

(02 " 02) 

(03 “ e 3) 


c ltQ 

°lt x 

°Z-t2 



In applications of the method of Iteration, It will 
be found that in many cases various approximations may be 
emnloyed quite advantageously. Thus, for the first and 
second iterations, the equivalent geometric angle-of- 
attack distribution can often be approximated by uniform 
and linear distributions, respectively, so that the 
c distribution may be obtained directly from the data 

of reference 7> which gives results for a wide range of 
taper ratios including the low aspect ratios commonly 
employed on tails. In other cases, an approximation to the 
cjt distribution can be obtained rapidly by the method 
given In reference 8. 

The fact that the method of iteration is based on a 
procedure In which the twist obtained from each iteration 
is used to Initiate the torque and the twist of the 
succeeding Iteration permits. In many cases, a rapid 
estimation of the twist distributions for the higher- 
order iterations. Thus, Inasmuch as the twist for a given 
torque distribution Is directly proportional to the 
magnitude of the torque, it follows that the propor- 
tionality of the twists obtained at a section In succeeding 
iterations will depend on the similarity for the two 
Iterations of the distributions of torque. (See equa- 
tions (11) and (12).) Because the shapes of the torque 
and hinge-moment distributions t(y) and h(y) 

(equations (7) and (8)) for a particular tail are 
usually not very sensitive to the spanwise variations 
of 0 and j2f, the shapes of the twist distributions 
tend to resemble the corresponding twist distributions 
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that initiated them. In these cases, therefore, the 
twist distributions for a higher-order iteration may be 
estimated from a knowledge of the values obtained in the 
preceding iteration by taking the ratio of the twists 
for two consecutive Iterations at any suitable reference 
station. The twist 0 at any station for the iteration 
of order n is then given by 

e n (y) = 

°n-2 

reference station 

A similar procedure may be followed to estimate the 
increments in # and cj t for the higher-order 

iterations . 

The foregoing iteration procedure leads to series of 
the form 

C U = C U 0 + c Ui + C U 2 + C1 tj + • • • 

9 = 0 q + 0 ^ + 0 2 + + . . . 

0 = flo + 01 + 02 + ^3 + • • • 

Prom the formulas derived in the preceding section, the 
quantities ®1* and $1 will be noted to contain 

the dynamic pressure q as a factor; cit 2 * and 0 2 , 

which are dependent on the corresponding values obtained 
in the preceding iteration, contain q2; and so on. The 
lift coefficient and twist at a section may each be 
represented, therefore, by a power series in q. The 
coefficients of these series, which depend on the various 
aerodynamic, geometric, and structural parameters, in 
general vary with speed because of modifications in the 
aerodynamic characteristics of the tail introduced 
principally through the effects of compressibility. 

By the application of these iteration procedures, 
the elevator contribution to the pitching moment about 
the airplane center of gravity is obtained from 
equation (1) as follows: 
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nbt/2 

-H4 7 (° ltd + 0 Hi + a H2 + * * 

Art/2 


) otdy + (Tq +.Ti t. T 2 . • •) 


J13) 


where the numerical subscripts refer to the order of the 
twist iteration and Tn = 2[t n (0) + h n (0)^J. (See 
equations (7) and (8).) Equation (13) may be written 

M = Mq + Mi + M2 + . . . (lij.) 

The elevator reversal speed is obtained from the 
value of the dynamic pressure q that makes the right- 
hand side of equation ( 1I4. ) equal to zero. 

In a similar manner, the elevator hinge moment for 
the flexible tail may be obtained as 

H = Hq + Ei + K 2 + . . . (15) 


where H n = 2h n (0). 

An iteration procedure similar to that described for 
the elevator effectiveness may be followed to determine 
the effect of angle of attack of the tall and expressions 
for M and H similar in form to equations (lij.) 
and (15) will be obtained. 

The tendency noted previously with regard to the 
similarity In the respective distributions for 9, 
and c for the higher-order iterations will lead In 

many cases to considerable simplification in the 
iteration procedure for determining the effect of tail 
flexibility on the' pitching moment M and the hinge 
moment H. The increments in M and H for the higher- 
order Iterations can therefore be obtained in these 
cases by means of the following relationships: 


*n = 


M_ -t 2 

_ **-n-l 


M n -2 


II 


H n = 


n-1 


H., 


n-2 


(16) 


(17) 
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The forces and resulting twists on the tail structure 
are directly proportional to the tail angle of attack and 
the elevator deflection corresponding to the values that 
would be obtained in an assumed rigid structure or 
to at R and 6 R , respectively. (See table III.) The 

differentiation of equations (llj.) and ( 15 ) with respect 
to at R and 6 R and conversion to the nondimensional 

form gives, therefore, for the flexible tail the parameters 

and ^£EL. 
de R* 6a tR d5 R 6a tR 


Corrections for Compressibility 

Inasmuch as the aerodynamic charterl sties of the 
tail are affected to an important extent by compressi- 
bility, the effects of compressibility must be considered 
In predictions for the control characteristics of the 
tali. In the absence of experimental data, the following 
corrections for compressibility, based on the theory of 
small perturbations, which is discussed In more detail 
In reference 9 * are summarized for the parameters 
Involved in the present analysis. These corrections may 
be applied at speeds below that at which the critical 
compressibility effects occur or up to a Mach number of 
approximately 0.60 In conventional airplanes. 

The span- load distribution In a compressible flow 
should be computed on the basis of a fictitious aspect 
ratio equal to the true aspect ratio, reduced by the 

factor \/l - M^, and the resulting values of cj^. 
obtained for this re duced f ictitious aspect ratio should 

be multiplied by l/tyl - M^. Thus, If the primes denote 
values obtained for the fictitious airplane, 

A t ' -Afc\/l - M 2 

and 

C H' 

c l* = 

u 


Vl - m2 
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The values for the parameters 




as obtained from low- speed data should he multiplied by 


the factor 


l/t - M2. 


The slope of the lift-coefficient curve in three 
dimensional incompressible flow is corrected for 
compressibility by multiplying it by the factor 


EA t + 2 

B = 

E'At' + 2 

■ 

in which the symbols E and E' represent the potential- 
flow correction for chord effect in incompressible and 
compressible flow, respectively. This correction applies 
specifically to elliptical plan forms but is approximately 
correct for other plan forms. In incompressible flow, 

E is the ratio of the semiperimeter of the ellipse to 
the span of the airfoil, as indicated in reference 10. In 
compressible flow, the ratio E' is that for the 
fictitious elliptical tail of span bt and aspect 
ratio At’. 

The derivative de/dcfy for compressible flow is 
computed on the basis of a fictitious tail length equal 

to the true tail length increased by the factor l/Vl - m 2, 
and the fictitious aspect ratios for the wing and tail 
equal to the true aspect ratios reduced by the factor 

VTTp. 


APPLICATION OP METHOD 
Data for Calculations 


Calculations' were made by the foregoing procedure 
of iterations for the effect of tail flexibility on the 
longitudinal control characteristics for two modern 
fighter airplanes designated airplanes A and B in 
order to illustrate the method and to obtain quantitative 
results for some typical cases. The computations were 
made, for both airplanes, of the tall effectiveness, the 
hinge-moment characteristics., - and the control-force 
gradients required- in recovery from ..dives at sea level 
and at an altitude of 30,000 feet. 


is 
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Figures 1 and 2 show the plan forms and dimensions 
of the horizontal tails for airplanes A -ana B, 
respectively. These figures also give the location of 
the* flexural axis as determined from stabilizer torsional- 
rigidity tests made in connection with the present 
investigation. The torsional -rigidity tests for air- 
plane A were made by the Langley Flight Research Division 
and for airplane B by the Langley Aircraft Loads Division. 
Both the stabilizer and the elevator for airplane A are 
metal covered, whereas airplane B has a metal -covered 
stabilizer and a fabric-covered elevator. 


The aerodynamic parameters for the two airplanes 
were based on low-speed data corrected for compressibility 
effects, essentially as described orevioucly. The basic 
data employed in tne calculations, the source from which 
these data were obbalned, and the compress? bill ty 
corrections applied are given in tables I and II. Average 
values along the soan were assumed for the puroi^eters 

; * 

H 


f d°h\ 

V> 3 R/ c 


and 



which were obtained on the 


basis of estimates for and 6Ch/d G t R from 

low- speed flight data. In the comoutatlons, the aero- 
dynamic centers of the tail sections were assumed to be at 
the quarter-chord ooints of the sections. 


The tail-stiffness data for the calculations were 
obtained from flexibility tests made on the stabilizer 
and elevator of the full-size airplanes. In order to 
clarify the relationship of the f lexibility-test results 
to actual flight conditions, the procedure for the 
determination of the stiffness data is described herein 
in some detail. The stabilizer tests were made by 
applying a concentrated torsional couple at a section 
near one tip of the stabilizer and measuring the 
torsional deflections at several stations along the span 
with reference to a station on the unloaded half of the 
stabilizer. The 6levator- flexibility tests were made by 
loading' bags containing lead shot or sand on one-half of 
the elevator along a line one -third of the chord behind 
the hinge with the elevator locked in position. The 
spanwise loading on the elevator surface corresponded 
approximately to a uniform distribution. The deflections 
of the elevator on the loaded side were measured at 
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several stations with respect to a reference station 
taken on the unloaded half of the elevator. 

Tests wer*e also performed- on one* rib of the- - ^ 
elevator of airplane B at a station about lj ..5 feet from 
the fuselage center line to obtain the effect on the 
distortion along the chord of a chordwise loading that 
simulated a triangular distribution more closely than 
the one just described (load concentrated at one- third 
of chord behind hinge). In these tests, measurements were 
taken of the deflections at very small intervals along 
the chord (5 dial gages for an 11 -Inch chord). The 
results Indicated that with both types of loading the 
distortions along the chord were equal and that the 
deflections along the chord .followed a straight line. 

It was assumed, therefore, that the measured angular 
deflection due to the elevator flexibility could be 
considered as an equivalent change in elevator deflection 
with no change in the camber of the elevator surface . 

The tail- stiffness data for airplane A are shown in 
figure 3 * The results for the stabilizer given in this 
figure are based on a concentrated torque of 833 foot- 
pounds applied to the right half of the stabilizer at a 
station 6.5O feet from the fuselage center line, and the 
data for the elevator are based on a total hinge moment 
of 83.3 foot-oounds distributed on the right half of the 
elevator in the manner described previously. 

The tail-stiffness data for airplane B are shown in 
figure Ij.. The results for the stabilizer given in this 
figure are based on a concentrated torque of 50U foot- 
pounds applied to the left half of the stabilizer at a 
station 5.92 feet from the fuselage center line, and the 
elevator data in figure ij. are based on a total hinge 
moment of 60 foot-pounds distributed on the right half 
of the elevator in the manner described previously. 


Procedure for Calculations 

With the aid of the foregoing data, the spanwise 
distributions for c 0 , and $ were determined for 

several Iterations. The oj t distributions were obtained 

by the usual methods based on lifting-line theory. In 
the computations, the stabilizers for the two airplanes 
were assumed to act in torsion similarly to tubes so that 


20 


NACA ACR No. L5B01 


the distributions of stabilizer twist resulting from the 
aerodynamic forces were calculated by means of 
equations (7) and (11) by use of the stabilizer torsional- 
rigidity coefficients shown in figures 3 and. 1;. Because 
the results of the elevator-flexibility tests for both 
airplanes indicated the probability that the static loads 
on the elevator did not act in pure torsion (see figs. 3 
and I4. for d^f/dH near root and tip sections), it was 
believed Dractical in the present investigation to modify 
the method described previously for determining $ . The 
twist distributions due to elevator flexibility were 
obtained, therefore, by multiplying the total hinge 
moment acting on each half of the elevator by the rigidity 

factors ^(y) shown in figures 3 and 1]-. This method 
dH 

for determining the elevator twist is strictly correct 
only if the loading on the elevator surface in the static 
tests simulated the loading in flight. For the present 
investigation, however, the error from this source is not 
expected to be important. Some computations with 
different assumed elevator flexibilities, which are 
discussed in the section entitled "Results and Discussion", 
indicate that the calculated results are not sensitive 
to reasonable variations in the elevator flexibility. 

The increments for c 0, and $ that were 

obtained in the various iterations were used to compute 
the pitching moment about the airplane center of 
gravity M by means of equations (13) and (lij.) and to 
compute the elevator hinge moment H by means of the 
corresponding equation (15)* The results for M and H 
were then converted into the nondimensicnal form as the 
derivatives dCm/dat R , £Ch/<bat R , and dCjh/de R . 

The details of the computations for determining 
and dCh/dSR for airplane B for a Maoh number 
of 0.60 at sea level are shown in table III. 

The change in control force per unit change in 
normal acceleration in recovery from a dive was computed 
by means of the following formula: 

K e qc Q 2 b e (18) 
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where ABr and Aat R refer to the changes In Br 
a'nd cL-t R per unit change 'In normal acceleration in- terms 
of g and where K e is the- elevator gearing ratio. In 
equation ( 18 ) 


A Or 


" 

— 

n de 


“ 

w 

dC]i] R ^ d 

Br\ “ do,, 

/dC m 

d c m R \ 

IV s 

^ 6Br 

dBj^ 

=4 ' a 

\ 4 °tR 

4a %/_ 


+ 28.6i t pg 


©a 


(19) 


and 


Aat T 


= - 


) + 28. 6 it Ps 

dOwy 


( 20 ) 


In equation (19) the two terms on the right-hand side 
enclosed In the brackets represent the part of ABr 

required to trim the airolane, and the third term 
represents the part of ABr required to balance the 
effects of rotation of the tail during the steady ohase 
of the pitching motion. In these equations, dC^/dOR, 

dCh/da tR , d>Cm/d6R, dCm/datj*, and are the 

values for the flexible tail obtained by the iteration 
procedure. If values for these parameters for the 
assumed rigid tail are used In equations (18), (19)* 
and (20), the control-force gradient for the rigid 
tail Pno is obtained. The values for the derive- 

<& 

flight results at an indicated airspeed of approxi- 

/dBjA 

mately 200 miles per hour. A value of equal 


tl ve 


for airplanes A and B were based on 


a 
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to -3.20 was obtained for airplane A based on a center- 
of -gravity location of 28 percent' of ‘the mean aerodynamic 
chord; whereas a value of - 3 .26 was obtained for air- 
plane B based on a center-of -gravity location of 
29.5 percent of the mean aerodynamic chord. The effect 
on Fn and F nR of movements of the airplane center of 

gravity was investigated by assuming different values 
for « . 

Wr 


The variation of 


de 
do w 


with speed as determined by 


means of the theoretical compressibility corrections 
noted previously, in conjunction with the design charts 
of reference 11, indicated a negligible change in this 
parameter up to a Mach number of 0.60. It was therefore 
considered sufficiently accurate in the present compu- 
tations to assume constant values for dc/do^. The 

values for , however, were corrected for compressi- 

'r 


bility effects by multi dying the low- speed value by the 
factor l/B corresponding to an average between the wing 
and tail. 


RESULTS AND DISCUSSION 


The results of the calculations are presented In 
table III and In figures 5 to 9* Table III shows the 
results obtained for the various Iterations In 
determining dCm/dbR and dCh/^bR for airplane B at a 

Mach number of 0.60 at sea level. Figure 5 shows the 
spanwise variation of tail angle of attack at and 
elevator deflection 6 resulting from an application of. 
the elevator control equivalent to unit deflection for 
the assumed rigid tall, as obtained from table III. 
Figures 6 to 9 shew the effects of horizontal-tail 
flexibility on the longitudinal control characteristics 
for airplanes A and B for a range of true airspeeds 
from 0 to 550 miles per hour at sea level and from 0 
to i +90 miles per hour at an altitude of 30>000 feet. 

This range of true airspeed corresponds to a Mach number 
range from 0 to 0 . 7 2 .. The speed for each altitude 
corresponding to a Mach number of 0.60, which represents 


NA0A ACR JNO. L5B01 


23 


the limit for which the theoretical compressibility 
corrections employed in the present computations are 
believed to be reliable, is indicated on figures 6 to. 
The results for the Mach numbers higher than 0.60 are 
included in the figures in order to give an indication 
of the trend of the flexibility effects. The results of 
the computations for a Mach number of 0.60 at sea level 
are summarized for both airplanes in table IV. 


Table III indicates that- the convergence of the 
iteration procedure is very rapid. This convergence, as 
indicated for dCm/dbR and dCh/d6 R for airplane 5 , 

is typical for the other parameters in the flexible tail 
for both airplanes. Thus, if the contribution obtained 
for each successive iteration is expressed as . a ratio of 
that obtained for the zerotli-order twist iteration. 


ac mR /66 R 


1 - 0.301 +0.0727 - 0.0167 + 0.00381+ = 0.759] 


>( 21 ) 


aCj^/d 6 R 
ach R /ab R 


1 - 0.21+1 +0.0557 -0.0131+0.00306 = 0.805 


The subsequent comparison illustrates the number of 
twist iterations required by the regular procedure of 
iteration In order to determine the longitudinal control 
characteristics for a flexible tail. The results obtained 
from table III as given by equations ( 21 ), which utilized 
three regular twist Iterations, will be compared with 
results obtained by the use of one and two regular twist 
iterations, respectively. In conjunction with the 
relationships given by equations (l6) and (17) for 
estimating M n and H n for the twist Iterations of 

higher order than one and two, respectively. Thus, by 
use of one regular twist Iteration, 


dC m / 65 R 

^Cxp-r/^Br 


dCh/dCR 

dChR/dbR 


1 - 0.301 +0.0905 - 0.0272 + 0.00815 = O.77I 


r(22) 


1-0.2^1 + 0.0581 -0.01^0 +0.00338 = 0.806 


2k 
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By use of two regular twist iterations. 




= 1 - 0.301 + 0.0727 - 0.0176 + 0.001+25 = O.759 


6C h A6 R 

= 1 - 0.21+1 + 0.0557 - 0.0129 + 0.00300 = 0.805 

dChp/o Bp 


(23) 


The comparison of the results shown in equations (22) 
with those given by equations (21) indicates that the use 
of one regular twist iteration in conjunction with the 
simple relationships given by equations (lb) and (17) is 
sufficient to determine the effect of tail flexibility 
to an accuracy of the order of 1 percent. 


Figures 6 and 7 show for airplanes A and B, 
respectively, the ratio of the tail effectiveness and 
hinge-moment parameters as obtained in the actual flexible 
tail to those obtained for the assumed rigid tail. These 
figures indicate, for the complete range of airspeeds for 
both airplanes, that the parameters t C^/6 Sp, dCm/da-fc^, 

dCh/£BR» and dCh/^ a tp are reduced numerically because 


of the tail flexibility and that the parameters 

dc h /da R 


/d 6 R \ 


and 


are increased numerica7i.ly because of this 


factor. The numerical reduction in caused by 


tail flexibility is due to the fact that the center of 
pressure of the lift resulting from the elevator 
deflection is behind the flexural axis (see figs. 1 and 2 
for flexural-axis locations) and the resulting torsional 
moment twists the stabilizer in a manner that reduces 
the tail lift. The numerical reduction in dCm/dBp is 

also due to the negative value of dCh/d5 R , which causes 

the elevator to twist and thus to reduce the elevator 
deflection. 
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The numerical reduction in dCm/dat R due to tail 

flexibility re suited, be cause the location of the flexural 
axis of the stabilizer is ahead of "the aerodynamic center 
and because thei value of dCh/dat R i a negative. The 

respective numerical reduction in the values for dCh/dCR 
and dCh/&at R due to the tail flexibility resulted 

principally from the fact that each of these parameters 
for the rigid tail is negative and' the elevator twist 
therefore numerically reduces the hinge moment in each 
case. The forward Dosition of the stabilizer flexural 
axis relative to the center of pressure of the lift 
contributed by the elevator tended, however, to increase 
numerically the value of due to the stabilizer 

twist (fi is increased by -0); similarly, the location 
of the flexural axis of the stabilizer ahead of its 
aerodynamic center tended to increase numerically the 
value for dCh/dat R . 


Figures 6 and 7 indicate that, in general, the 
effects of tail flexibility vary with speed and altitude 
approximately as the dynamic pressure - modified, of 
course, by the relative compressibility effects. This 
variation with speed and altitude results from the raoid 
convergence of the power series in q, which causes the 
terms in q of higher order than unity to be compara- 
tively small. In some cases, however, at very high 

/ dCh/datp 

soeeds [see figs. 7(b) and for — — and Fn/Fn R j 

\ i>CbR/W R 

respectively), the effects of the terms in q of higher 
power than unity become comparatively significant. 


Computations were made to estimate the effect on the 
parameters shown in figures 6 and 7 o f increasing the 
elevator stiffness at each section by 12.5 percent of the ■ 
average elevator stiffness. The results of these 
computations indicated that, for a Mach number of 0.60 
at sea level, the ratios of the parameters dCm/dOR, 

dCh/dBR, and dCh/dat R to the corresponding ratios for 

the assumed rigid tail would be increased in the order 
of 2.5 percent as compared with those shown in figures 6 
and 7, and the corresponding ratio for dC]n/dat R would 

be increased by less than 1 percent; whereas, at 
30,000 feet for the same Mach number, the effect of the 
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Increased elevator stiffness would be about O.ij.0 of the 
corresponding foregoing effects indicated at sea level. 

It can be noted from figures 6 and 7 that, provided 
critical compressibility effects do not appear, elevator 
reversal for both airplanes A and B does not. occur up to 
a speed corresponding to a Kach number of 0 . 72 . 


Figures 8 and 9 present a comparison of the control- 
force gradients in recovery from dives as obtained for 
the actual flexible tail and assumed rigid tail. It 
should be noted in these figures that the required motions 
of the elevator control stick per unit g are not 
necessarily equal for the flexible and assumed rigid 
tails. Figure 8 gives the results for airplane A at sea 
level and at an altitude of 30,000 feet. This figure 


shov/s the variation with airspeed of 

/ d6 R\ 

ratio Pn/Pn R for values of f — — J 


F n and the 
in Incompressible 


flow of - 3.2 and -1.6o. These values of - 3.2 and - 1.60 
correspond, respectively, to center- of-gravitv locations 
at 28 percent and approximately 31 percent of the mean 
aerodynamic chord. Figure 8 shows that flexibility of 
the tail increases the control-force gradient and that 
this increase for a Mach number of 0.60 amounts to 


12 percent at sea level and 3*5 p^conl at 30,000 feet 
altitude . This figure also shows that a rearward movement 
of the center of gravity of approximately 3 percent of 


the mean aerodynamic chord causes a small reduction in 
the ratio F-Zlv, . The results for the airplane B at 

R /aei 

sea level and at altitude for values of — \ in 


incompressible flow of - 3.26 and - 6.00 are presented in 
figure 9 . These values of -3*26 and -6,00 correspond, 
respectively, to conter-of-gravity locations at 29*5 percent 
and approx imately 25 percent of the mean aerodynamic 
chord. The figure shows, for airplane B for a range of 
airspeeds at the altitudes considered, a small Increase 
in the control-force gradient due to tail flexibility, or 
approximately one -half of that Indicated in figure 8 for 
airplane A. Figure 9 also shows that a forward movement 
of the center of gravity of approximately l \., 5 percent of 


the mean aerodynamic chord causes a small Increase in the 
retio F n /F nR * 
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An examination of equations (l8), ( 19 ), and (20) 
Indicates th£t ' the control -force gradient ...in_ a ..dive 
recovery may be Influenced to an Important extent by- the 

^Ch/o8R 

aerodynamic parameters '* ^Ch/^ a tRj and dCm/oat^. 


The results of the present analysis show for both air- 
planes A and B that the first two of these parameters are 
affected by tall flexibility in a manner to Increase F n ; 
whereas dCm/dat^ is affected by this factor in a manner 


to reduce F n . As noted previously, the numerical reduction 
in dCm/batR obtained in the present computations for 


airplanes A and B is caused principally by the location 
of the flexural axis of the stabilizer ahead of its 
aerodynamic center and by the negative value 
of dCh/datR* In order to obtain an Indication of the 

importance of the change in dCm/dat R due to tail 


flexibility for the control-force gradient in a dive 
recovery, commutations were made for the two airplanes In 

which it was assumed that which is roughly 

da tB <ba tR 


equivalent in the present case to a rearward movement of 
the flexural axis back to the aerodynamic center. These 
computations indicated, for a Mach number of 0.60 at sea 
level, that in the case of airplane A the ratio Fn/F nR 

would be increased from 1.12 to 1.26, and In the case of 
airplane B this ratio would be increased from 1.0 5 
to 1.075* On the basis of the present analysis it appears, 
therefore, that the location of the flexural axis of the 
stabilizer too far behind the aerodynamic center of the 
tail, could cause excessive control forces in a dive 
recovery at high speeds. 


CONCLUSIONS 


An Iteration method for determining the effect of 
tail flexibility on the longitudinal control charac- 
teristics of airplanes was applied to two modern fighter 
airplanes and was found to provide a practical procedure 
for the determination of these effects. 
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The re suit 3 of calculations to determine the effect 
of bail flexibility on the longitudinal control charac- 
teristics for two fighter airplanes Indicate that the 
longitudinal control characteristics are affected to a 
significant extent at high speeds by this factor. The 
following conclusions apply to re3ult3 for these airplanes 
at sooed3 below that at which critical compressibility 
effects occur: 

1. Tho magnitude of the tall-f?.e;cibllity effects. 

In general, varied approximately as the dynamic pressure - 
modified, of course, by the relative compressibility 
effects. In some cases at very high speeds, however, 
the effects of. the terms containing the dynamic pressure 
of powers greater than unity becune comparatively 
significant. 

i 

2. Tail flexibility was found to reduce significantly 
the rates of change of pitching moment and hinge moment 
with elevator deflection and tail angle of attaok. 

3 * The control-force gradients In a dive recovery 
were Increased because of tail flexibility. 

ij.. Rearward movements of the airplane center of 
gravity tended to decrease the effects of the tail 
flexibility on the control-force gradient; whereas 
forward movements of the airplane ■ centor of gravity 
tended to increase the magnitude of these effects. 

5. The location of the flsxural axis of the 
stabilizer relative to the aerodynamic center of th9 
tail is an Important dasign consideration with regard 
to the magnitude of the ■ tail-f loxibility effects. The 
location of the flexural axis of the stabilizer too far 
behind the aerodynamic canter could cause excessive 
control forces In a dive recovery at high speeds. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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tab lx i 

DATA FOR CALCULATIONS - PHYSICAL AND GEOMETRIC CHARACTERISTICS 
[pat a rural ah* d by natmf acturer] 


Airplane 

Weight, 

W 

(lb) 

Wing 

area, 

S 

(»q ft) 

m 

Mean- 

aero- 
dynamic 
chord of 
wing. c. 
(ft) 

Tail 

area, 

(sq ft) 

Elevator 

•pan* 

b e 

<•« ft) 

Root mean 
chord or 
elevator, 

(ft) 

Tail 

length, 

(fti 

(radian/rt) 

A 

12,000 

300 

5.55 

7-28 

55 

16 

Wm 

21.4 

0.66 

B 

7*660 

2}6 

5.815 

6. 6J, 

lti.i 

ij.2 

Efl 

15.5 

-57 


TABLE II 

DATA FOR CALCULATIONS - AERODYNAMIC PARAMETERS 


Parameter 

Value in 
Incompressible 
flow 

Source of data 

Correction for 
compressibility 

Airplane A 

(6o b /66r) C i ^ 

*-0.0086 

Reference 12 

Multiply by r -k J 

“0 

0.095 

. Assumed 

Do. 


-0.00686 

Unpubli.h.d d.t. b.a.d on 
dC h /6otp * -0.00218 wnd 
dCh/dO R * -0.00801* 

Do. 


-0.0552 

None 

( ) 

«-0.66 

Reference 12 

None 

d c/da* 

0.50 

Reference 11 

Assumed constant 


0.077 

Reference ^ 

Multiply by B 

(«r/<» c l) r 

-J.2 

Based on unpublished date 
for c.g. at 28 percent M.A.C. 

Multiply by 1/B; 
average for wing and tall 

Airplane B 

gggppai 

8-0.0091 

Reference 12 

■BSXS&3II 

*0 

0.095 

Aeeumed 

Do • 

(dc h /d8 R ) oit 

-0.00605 

Estimated from unpublished, 
flight data based on 
dC h /dat R = -0.000511 and 
dG h /d8 R = -0.00655 

Do. 

(6Ch/6o lt ) SR 

-0.00825 

None 

(da tR /60 R ) 

l t 

8-0.59 

Reference 12 

None 

dt / d ‘‘w 

0.50 

Reference 11 

Assumed constant 

a 

0.077 

Reference 7 

Multiply by B 

(dO R /dC L ) R 

-3.26 

Estimated from unpublished 
flight data for c.g. at 
£9.5 percant K.A.C. 

Multiply by 1/B; 
average for wing end tall 


“Value given la for a section at 4.5 ft. from fuselage center line; appropriate values were used for other 
sections . 

^Average constant values were used. 
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TABLE III.- COMPUTATI OKS FOR 57S AND OF AIRPLANE B AT A MACH NUMBER OF 0.60 AT SEA LEVEL 

R i5 R 

jcompreael bl 11 ty effecta have been accounted for in all aerodynamic parameters; 


q - 53U. 6 


/dC m \ /dC>A /dc h \ 1 

: (do?) = -°- 01 ^ dS? = (5jf-) = -0.00325 

v v «/= lt V Ol ty 0R J 


If 



Distance 

from 

center 

line 

(ft) 

0 

a*. 

(deg) 

c l 

(tase 5 

on 

eq. (2) 

e 

(n ? . 

=t 2 

dT 

d U 

(eq.(3)) 

4 

db 

dr\ 

(eq. (5) ) 

t 

(eq. (7) ) 

°TR 3 

(f 4?' 

0 

(eq.(ll)) 

h(0) 

(eq.(S)) 

dH 

(fig-4) 

*< 

=t 

nr 

From table I , 

i t = 15.5 ft, s = 236 sq ft, c w = 6.64 ft 
C 9 = 1.01 ft, b 8 - 13.2 ft 

Zeroth order of twist Iteration 

* -2l t q l r 2 \ P r -| 

6 

nr 

C1 t 0 

6 S 

e 

c t 

c t 2 

d ' f o/ d n 

c 2 
c e 

dhQ/dr) 

§ 

c tr 3 

e l 

h 0 (0) 

dg 

dH 

01 

5 r 

c t 

5 r 


»R 

= -15,660 ft-lb 

S'T" -59.03 rt-ib 

0-575 

1.25 

2.50 

5-75 

5.00 

5.85 

6.60 

0.575 

1.25 

2.50 

5-75 

5.00 

lil 

0.575 
1.25 
2.50 
5-75 
5.00 
5-63 
6. 60 

r 

1 


0.0425 

.0440 

.0455 

x 

.0342 

0 



0.026 

.00S 

-.010 

-.036 

-.068 

-.116 



16.6 

14.5 

12.0 

M 

6.55 

-91.36 

-Ml 

-67.74 

-57.31 

-52.17 

I.89 

1.6k 

1-37 

1.00 

0 

-5.812 

-4.241. 

-3.251 

-2.598 

i231.2 

-147,6 
-71+-17 
-31.01 
; 0 

00 

21200 

10200 

5560 

2560 

1800 

0 

-.00791 

-.03229 

-.06362 

-.096^ 

-.11536 

-.12329 



-29.51 

0.0063k 
. 00700 
.00835 
.OHIO 

.01519 

.01U50 

-0.1871 

-.2066 

-.2464 

--■m 

-.4279 

4.08 

3.81 

3-47 

5-12 

2-75 

2.52 

0 

First order of twist lteretlon 

%-'-^r ([ o 2 v td T* IS' 01 + h i (oi l 

= 4711 n-16 

H. 2h,(0) 

Rr= Sr = 14.25 n-ib 

ft - a i 

L 

°R 

° l t 1 

e 

ct 

=t 2 

dT-i/dn 

*. 2 

dh^/dri 

°R 

k 

5, 

c tr 5 

m 

h](0) 

d 0 

dH 

£2 

5 r 

C t 

°p. 

«R 

5 r 

6= 

-0.1871 

-.1937 

-.2141 

-.2640 

-0.1671 
-.2066 
- .21+6^ 

'-Ml 

-.4279 

-a00996 

-.01074 

-.01277 

-.01526 

-.01599 

-.01547 

0.026 

.008 

-.010 

-.036 

-.068 

-.116 

16.6 

14-5 

12.0 

9-73 

f.56 

6.35 

0 

16.66 
16.92 
16.46 
18. 52 
20.58 
17.38 

HI 

1.37 

1.00 

:W 

0 

nil 

1.258 

1.143 

!.i43 

.829 

100.72 

86.88 

67.16 

46.76 

23-21 

7.89 

00 

21200 

10200 

5560 

296° 

1800 

0 

. 00157 
.00886 
.01841 
.02907 

■m 

7-116 

O.OO65I* 

.00700 

.00635 

.01110 

.01519 

.01450 

0.0452 

■0497 

.0505 

.0789 

.1080 

.1032 

4.08 

3-81 

3-47 

3.12 

2.75 

2.52 

0 











Second order of twist Iteration 

/h \ 

I 0 'HH + 4^ ( ° )th2,0) ] 

= -1139 ft-lb 
Ho 2ho(0) 

= — Sa-L = -5.29 ft-lb 
°R ’r 

02 - 3 2 

0 2 

6 R 

IT 

e 

c t 


dTg/dn 

°e 2 

; 

dh^d*^ 

t^ 

6 r 

c tr 3 

5i 

S R 

h 2 (0) 

dg 

dH 

tl 

»R 

°t 

°R 


5 r, 

^R 

5 r 

0.04516 

.01,776 

.05059 

.05762 

A 

■0595 

.0789 

.1080 

.1052 

0.002443 

.00261k 

.003106 

.003719 

.OO39I4 

.OO3285 

0 

0.026 

.008 

-.010 

-.036 

-.068 

-.116 

..... 

16.6 

14.5 

12.0 

V.ll 

6.35 

L! 

-4.007 
-4.055 
-3- 9°i 
-4.H4 
-4.714 
-3-959 

I.89 

1.64 

1-37 

1.00 

:L 65 

0 

-0.3652 

-M 

-.2497 

-.2593 

-.1813 

-23. 16 
-19-88 
-15.30 
-10.60 

0 

00 

21200 

10200 

5560 

2560 

1800 

- ■ ■ ..j 

0 

-.OOOkl 

-.00161 

-.00k06 

-.00622 

-.00775 

-.00B21 

-1.61*3 

0.0063k 

.00700 

.00835 

.OHIO 

■01519 

.01450 

-0.0101+2 

-.01150 

-.01372 

-.0182k 

-.02496 

-.02382 



4.08 

3.81 

5.67 

5.12 

2.75 

2.52 

0 



NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


w 

to 


NACA ACR No. L5B01 











60m dCv, 

TABLE III.- COMPUTATIONS FOR T7~ AND — - OF AIRPLANE B AT A MACH NUMBER OP 0.60 AT SEA LEVEL - Concluded 

“°R OOR 

Third order of twist Iteration 


By procedures similar to those given previously 


Mi 

•P = 261.6 f 

6 R 

h 3 

x 2 - = 0.771 f 


771 ft-lb 


Fourth order of twist Iteration 


By means of equation (16) 


3 * - - 60 . 


08 ft-lb 


By means of equation (17) 


ft- 


1807 ft-lb 


If equations ll* and 15 are employed and the results expressed nondlmenslonally 


, (*0 + *1 * I *3 + \)/ a H = _ 0i0 _ 

d8 R qsc_ 


6y Sc » _ 8 yds R 
M 0 " d <V>®R 


= 0.759 


*Jh = (H 0 + H x + H 2 + H, ♦ H L )/S R 
d5 R <l~ e \ 


= -0.00660; 


°°h 

8C l/ d5 R 
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TABLE IV.- COMPARISON OP EFFECT OF HORIZONTAL-TAIL 
FLEXIBILITY ON LONGITUDINAL CONTROL CHARACTERISTICS 
FOR AIRPLANES A AND B AT A MACH NUMBER 
OF -0.60 AT-SEA- LEVEL 


Parameter ratio 

Airplane A 

Airplane B 

dC m /&6R 

0.79 

O .76 

dCnj/da tR 

dC mR yka tR 

■97 

• 98 

(6 6ly /4a tR ) to 

1.22 

1.30 

6Ch/6GR 

6Chjy'65 R 

• 91 

j 

.30 

dC h /da tR 

V* 1 ! 

• 9k 

.96 

6Ch/^5R 

6Cuj/^6 r 

d ° h R/ fe6 R 

^C mR ^b6 R 

1.15 

1.06 

*n_ 

Pn R 

1.12 

1.06 
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Distance from ta/I center ! me ? ft 
Figure 3.~ £ x peri mental data for flexibility ° 
hor/ziontal tail. Airplane A. Data from tests 
made by Langley Flight Fe search Dn/ision. 






Stabilizer tunst ? 6 } deg Elevator r/giditg ; d(p/dli ; deg/ ft- lb 
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Distance from tail center hne ? ft 

Figure 4.~ Experimental data for flexibility of 
horizontal tail . Airplane B . Data from tests 
made by Lang /eg Aircraft Loads Division. 
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Fig. 5 



D tsia nee. from fa/7 center hne^ ff 


F/gure S- Distribution of fail ang/e of affacf ond e levator 
deflection m a flexible, fail resulting From an appli- 
cation of l he elevator control equivalent: to unit 
deflection for the assumed rigid toil. A irpfane. 3j 
Mach numbegOEOat sea level OC tf ^0 0 . 


le of attack; oc^, o/e. 
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Altitude 

•U 


ft) 


— 

Sea level 
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— 

-- 3 OP 00 


*| Indicates limiting speed 
for u/h/ch calculated 
compressibility corrections 
re believed reliable 


True airspeed^ mpf? 

(d) Effectiveness parameters . 

F/pure 6 -Effect of horizontal -tail flex/bil/ti / 
or longitud '/ na I -control parameters. Airplane A. 


ig. 6a 
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0 /OO ZOO e 300 400 S00 600 

True airs peed ^ rn p h 


Id) E ffec it veness par a me te rs. 

F/pare 7- Effect of horizontal-tail flexibility 
on lonpi tud/na/- control parameters. f\ ir plane B. 
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NACA ACR No. L5B01 


^Pi S • 7b 
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.2 

(9 

O IOO 200 JOO 400 SOO 600 
True airspeed; m p d 

(b) H mge- moment parameters . 
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*1 Indicates limiting speed 
; for which calculated 
compressibility corrections 
are believed reliable- 
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Figure 7 - Concluded. 
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F n for flexible tail 
F nfl for rigid tail - 
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Fig. 9a, b 




&)Atsea /eve/. (b) At ofirtude of 30,000 feet. 

F/yar& ‘/.-Effect of horizontal - 1 a// flexibility on elevator 
control- force yrad/ents m recovery from d/ve.. 
Airplane B. 
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